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We have performed the numerical modehng of Ge(lll)-(2xl) surface electronic properties in 
vicinity of P donor impurity atom located near the surface. We have found a notable increase of 
surface LDOS around surface dopant near the bottom of empty surface states band tt*, which we 
called split state due to its limited spatial extent and energetic position inside the band gap. We 
show, that despite of well established bulk donor impurity energy level position at the very bottom of 
conduction band, surface donor impurity on Ge(lll)-(2xl) surface might produce energy level below 
Fermi energy, depending on impurity atom local environment. It was demonstrated, that impurity, 
located in subsurface atomic layers, is visible in STM experiment on Ge(lll)-(2xl) surface. The 
quasi- ID character of impurity image, observed in STM experiments, is confirmed by our computer 
simulations with a note that a few vr-bonded dimer rows may be affected by the presence of impurity 
atom. 

PACS numbers: 68.35.Dv, 68.37.Ef, 73.20.At, 73.20.Hb 



I. INTRODUCTION 

At present it is a common place that Ge(lll)-(2xl) 
surface consists of 7r-bonded zigzag chains. This was 
confirmed many times by different means, se^^H^ and ref- 
erences therein. Surprisingly, just a few publications are 
devoted to investigations of impurity atoms on (111) sur- 
face of elemental semiconductors^^^. And as one can 
see the interest to Si(lll)-(2xl) surface is renewed. But 
not to Ge(lll)-(2xl) surface. This is unexplainable, be- 
cause Ge is the main candidate for technology, allowing 
to overcome scaling limits of Si-based MOSFETJ^. The 
knowledge of local properties of Ge, especially caused by 
impurity atoms, is of vital importance. Besides, the sur- 
face and interface properties of Ge(lll) have great sig- 
nificance for Ge spintronics applications. It is known Ge 
has some advantages above Si^^. 

The STM method till now is the only physical method 
achieving atomic resolution in real space. But exper- 
imenters often suffers from the lack of some reference 
points provided by the theory. For example, the reliable 
STM image interpretation still remains the challenging 
task. There is no general approach which takes into ac- 
count all kind of physical processes responsible for STM 
image formation. Below we report on surface electronic 
structure investigation performed by ah initio computer 
simulations in the density functional framework, which 
is the first order estimation for STM/STS images, and 
could serve as a basis for further model improvements. 

We restrict the present investigation to the case of left 
(negative) only surface buckling (Fig.jljl)) as at present 
this matter is still controversial and is the subject of in- 
tensive invest igat ion Our research is in some sense 
similar to reported in^for Si(lll)-(2xl) surface. But we 
take into account all possible impurity positions in two 
surface bi-layers of Ge(lll)-(2xl) reconstruction, and be- 
sides, our analysis is not aimed on pure STM image simu- 



lation, rather on comprehensive analyses of local density 
of states. 



II. METHODS 

We have performed our DFT calculations in LDA ap- 
proximation as implemented in SIESTAI^ package. The 
use of strictly localized numerical atomic orbitals is neces- 
sary to be able to finish the modeling of large surface cell 
in reasonable time. The surface Ge(lll)-(2xl) super-cell 
consists of 7x21 cells of elementary 2x1 reconstruction, 
each 8 Ge atomic layers thick (total 2646 atoms). Vac- 
uum gap is chosen rather big - about 20 A. Ge dangling 
bonds at the slab bottom surface are terminated with H 
atoms to prevent surface states formation. The geome- 
try of the structure was fully relaxed, until atomic forces 
have became less then 0.003 eV/A. More details about 
calculations can be found elsewher^^. 

As we have reported earlier, the atomic structure of 
Ge(lll)-(2xl) surfa ce is strongly disturbed in vicinity of 
surface defect 4^^^. A few 7r-bonded rows around the 
defect are affected. That is why the geometry relaxation 
has been performed with the large super-cell to keep the 
internally periodic for DFT images of impurity well sep- 
arated. Although it is still the open question, if the de- 
fect's images separation is sufficiently large. 

At the last step of simulation the spatial distri- 
bution of Khon-Sham wave- functions and correspond- 
ing scalar field of surface electronic density of states 
LDOS{x^ eV) were calculated. Because of strictly lo- 
calized atomic orbitals, used in SIESTA, the special pro- 
cedure of wave- functions extrapolation into the vacuum 
has to be used (it is also implemented in SIESTA pack- 
age). 
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III. RESULTS 

A. Geometry and ground state properties 

In STM method the LDOS is measured above the sur- 
face, The tails of wave- functions are actuahy making the 
image. Thus in DFT calculations we are interested in the 
following quantity: 

LDOS{x, y, eV) - ^ \^{x, y)\'' 5{E - i?i)U=Con.t, 

where ^ are Khon-Sham eigenfunctions, 6 is finite width 
smearing function, Ei are Khon-Sham eigenvalues, and 
summing is evaluated at certain plane {z = Const), lo- 
cated a few angstroms above the surface. Here the broad- 
ening is the essential part of calculations, as we know 
from our experience that tunneling broadening in STM 
experiments on semiconductors typically amounts about 
100 meV. Broadening provides the degree of LDOS 
smoothing, necessary to resemble experimental tunnel- 
ing spectra. 

At the first stage of DFT calculations the equilib- 
rium geometry has to be established in the unit cell of 
Ge(lll)-(2xl) reconstruction. Afterwards the unit cell 
is enlarged to the desired extend, the defect is intro- 
duced and the structure is relaxed again. The final step 
is LDOS{x, eV) calculation. The results are sketched 
in Fig[l] 

The 1x1 and 2x1 surface irreducible Brillouin zones 
(IBZ) together with special points and directions are 
shown on pane (I). The right half of pane (I) contains the 
LDOS image correctly oriented with respect to special 
directions. The image corresponds to the whole 7x21 
super-cell of 2x1 surface unit cells of Ge(lll) surface, 
with P impurity atom placed at the position 1 (see be- 
low). We have to note that all figures are intensively 
cross-referenced, and the meaning of some notations on 
one figure can be better understood with another figure 
in mind. 

Pane Fig. [ijll) illustrates electronic structure of clean 
Ge(lll)-(2xiy surface. Two surface states (SS) bands, 
empty tt* and filled tt, can be seen in projected band 
gap. The widths of SS bands, derived from Fig. [ijll), 
equal to Att* = 1.24 eV and Att = 0.44 eV respectively. 
The surface band gap AEsbg is about 0.3 eV. As it 
should be expected, LDA approximation gives the band 
gap value which is much smaller than the experimental 
one. SS bands, as well as bulk bands, are also shown in 
Fig. [ijll) by small rectangles next to ordinate axis. This 
representation will be used on majority of figures below. 

The surface band structure is presented for the case 
of Ge(lll)-(2xl) surface with negative buckling. Our 
calculations predict that negative surface isomer is en- 
ergetically (by almost 11 meV per (2x1) unit cell) more 
favorable and, as we stated above, we will restrict present 
analyses to negative buckling only. 



The valence band (VB) top and empty SS bottom co- 
incide with Fermi energy in Fig. [ijll). The band gap is 
completely covered by empty SS tt* band. Specific for 
such band diagram was discussed in^. 




FIG. 1. (I) Sketches of left and right isomers of Ge(lll)- 
(2x1) surface. Two bonds, appropriate for definition, are 
marked by arrows. IBZ for Ge(lll) 1x1 and 2x1 surfaces 
together with special points and relevant directions denoted. 
Model surface slab with impurity atom positioned at site 1 
(see below) correctly oriented with respect to crystallographic 
directions. (II) Surface band diagram for Ge(lll)-(2xl) re- 
construction. Empty TT* and filled tt surface states bands 
can be seen in projected band gap. Their energetic posi- 
tion is schematically shown at ordinate axis by rectangles. 
To illustrate the relation between surface band structure and 
LDOS{eV) curve (III), the latter is shown at the left pane, tt* 
bottom and tt top are marked by arrows. Bulk bands energetic 
position is shown at the right side. (Ill) The LDOS{eV) 
curve and I(V) curve on logarithmic scale, averaged above 
the whole 8nm x 5nm surface slab shown on (I). 

Jumping ahead, on pane Fig. [ijlll) we display 
LDOS{eV) curve. The calculated I(V) dependence 
on a logarithmic scale is shown alongside. Surface 
LDOS{eV) curve is the result of averaging over the whole 
8 nm X 5 nm area on Figjljl). The gap right below the 
Fermi level is clearly observed on LDOS{eV) graph. It 
does not corresponds to bulk band gap. The closest re- 
semblance can be found with the surface band gap. 

Rectangles at absissa axis illustrate different band's 
energy position. It is necessary to state, that everywhere 
below, the bottom of CB is schematically shown on the 
figures for the case of Ge(lll)-(2xl) surface at room tem- 
perature. In this case the optical band gap is about 
0.5 eV^^. The DFT band gap in LDA approximation 
is non-physically small, less then 100 meV. To give clear 
impression on relation between surface band structure 
and LDOS, the latter is also depicted at the left side of 
surface band diagram Fig. [ijll) . 
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B. LDOS scalar field representation 

In further exposition we will focus mostly on LDOS 
properties and before we go to the main results we have 
to clarify the physical meaning of our data representation 
for LDOS. 



^-^ P at pos 1 1 . ^ 
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FIG. 2. (I) Sketch of LDOS{x,y,eV) scalar field. Rel- 
evant directions are shown. (II) Labels to identify donor 
atom position in two surface bi-layers of Ge(lll)-(2xl) re- 
construction. P atom is shown in position 1. (Ill) Quasi-3D 
and 2D representations of LDOS{x,y)\ev^o cross-section of 
LDOS{x, y, eV) scalar field at zero bias. Specific points and 
directions, referred to below, are denoted. 



Below we are speaking about cross-sectioning of 
LDOS{x, y, eV) scalar field (Fig. ^V)). The x and y di- 
rections correspond to [Oil] and [211] crystallographic 
directions. The two most relevant quantities are cross- 
section of scalar field LDOS{x^ eV) along (x, y) and 
(x, eV) planes - LDOS{x^ y) and LDOS{x^ eV) respec- 
tively. 

The LDOS is build up for different impurity atom po- 
sitions in two subsurface bi-layers of Ge(lll)-(2xl) re- 
construction. The definition for donor atom position's 
notation is presented in Fig. [2|ll). Worth to mention, 
that impurity atom does not occupy exactly the same 
lattice site as host atom it substitutes. 

In Fig. [2|lII) we show the cross-sections 
LDOS{x,y)\ev=o at fixed bias voltage for the case 
of P donor atom located at position 1 in surface bi-layer. 
These images roughly correspond to experimental STM 
images, as at small bias voltage there are not too many 
sharp LDOS features, contributing to the image. It 
can be seen from Fig. [2j that two 7r-bonded rows of 
surface reconstruction are influenced by impurity. In 
each row a protrusion can be observed. The spatial 
extent of impurity induced feature along the direction 
of TT-bonded dimer row ([Oil] direction) is at least 40 A. 
Note two distinguishable maxima on the protrusion. We 
will come back to this fact later. Arrows and dots on 
the image, as well as (a-a') and (b-b') lines, mark spatial 
points and directions, referred to on figures below. 

We have found notable increase of LDOS in vicinity of 
impurity atom at the bottom of empty SS band tt*. We 
will refer to it as to split state. Fig. [S] proves the valid- 
ity of this terminology. In the figure the LDOS{x^ y^ eV) 



field is shown by surfaces of equal value, colored by ap- 
plied tunneling bias voltage. The LDOS is drawn for two 
TT-bonded rows denoted as a-a and b-b' in Fig. jSJ The 
position of donor atom is depicted in the figure. To pre- 
vent the confusion caused by quasi-3D picture, we clarify 
that the impurity is located at the left side for 7r-bonded 
row b-b', and at the right side for 7r-bonded row a-a'. 

Let us point out some important facts about this 
LDOS representation. First of all, high values of LDOS 
are confined within areas bounded by surface of constant 
value. They are perfectly localized above up 7r-bonded 
rows. This is the reason, why only every second dimer 
row is imaged by STM. In between 7r-bonded rows LDOS 
is relatively low. Beside this, all round shaped vertical 
structures in 7r-bonded row are located above up-atom. 
Down-atom can be found in between them. Such spatial 
structure of LDOS is a consequence of collective 7r-bonds 
formation. STM can only image the up-row dimers, and 
therefore, only up-rows. Basically, with used approach 
TT-bonds can also be directly visualizecP^. 

The hybridization of atomic orbitals is clearly visible 
from Fig. [3j It is very strong in close proximity to surface 
defect. This causes the appearance of specific feature 
near the bottom of empty SS band (marked by arrows 
in Fig.|2|. The limited spatial and energy extents of this 
feature as well as its position in the band gap implies 
that this indeed is the split state. 



split state 



a_a. split state 




FIG. 3. LDOS{x,y, eV) field above rows a-a' and b-b' (see 
Fig. [2| shown by surfaces of constant value. Coloring corre- 
sponds to applied bias voltage. Split state in the band gap can 
be clearly observed above both dimer rows. Note the changes 
of structure at the right side of CB features and strong hy- 
bridization of atomic orbitals around impurity atom. 

We are working with microscopic picture, on the level 
of individual atoms. That is why we are able to see the 
connection between atomic orbital's hybridization and 
macroscopic band structure. Surface states itself appear 
due to atomic arrangement of the surface. Split state 
appear due to changes of this arrangement around the 
defect. They both have the same root - hybridization of 
atomic orbitals. One problem exists. It is really difficult 
to define the energetic position reference level, is it Fermi 
energy, or the SS tt* band bottom? To be accurate we 
will refer to Fermi level, i.e. split state is located at the 
Fermi level, and not near tt* bottom. 

Important for our result's understanding cross-sections 
LDOS(x^eV)\a-a'{b-b') fixed y coordinate are shown 
in Fig.pTI). They are taken along (a-a') and (b-b') planes 
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FIG. 4. (I) 2D LDOS{x,eV) distributions taken along 
a-a' and b-b' planes in Fig. [l] Impurity atom is at 1 site. 
Two rectangles at the left side correspond to optical and DFT 
conduction band positions. (II) Profiles of a-a' and b-b' 
panes along c-c' line, which are essentially the same as profiles 
of images in Fig. [l] along a-a' and b-b' lines. (Ill) Profiles of 
a-a' and b-b' images along d-d', e-e', f-f and g-g' lines, which 
are LDOS{eV) dependencies at points, shown in Fig. [l] by 
black dots and white arrows. For details see the text. 
Axis directions and images size are indicated on the figure. 



in Fig. [2j i.e. along 7r-bonded rows of Ge(lll)-(2xl) re- 
construction. Areas near the Fermi level, where split 
state resides, are zoomed in on the insets. The positions 
of Fermi level Ep, conduction band (CB) bottom, valence 
band (VB) top and empty tt* and filled tt surface states 
bands are indicated in Fig. [3j According to our DFT cal- 
culations the top of VB almost coincide with the bottom 
of empty surface states band tt* and Fermi level Fig. [l] 
an(i^. Even more, the tt* band can be partially filled 
at very high doping raticPI'. Our calculations also pre- 
dict for Ge(lll)-(2xl) surface with negative buckling the 
position of empty SS tt* band bottom about tens meV 
below the Fermi level. 

The proportions of LDOS{x^ eV) images are chosen on 
purpose in a way that is convenient for experimenters. 
Typically the number of point along spatial direction is 
less than the number of bias voltage points and tun- 
neling spectra image is elongated in vertical direction. 
Though the DFT band gap in LDA approximation is 
non-physically small for the sake of completeness it is 
also shown in Fig. |4jl). 

The distribution LDOSix^ eV) is reach of features. 
Its cross-section along x coordinate gives the LDOS{x) 
profile exactly in the same way as the cross-sectioning 
of LDOS{x^ y) (Fig. [T]) along x coordinate. The cross- 
sections of a-a' and b-b' panes along c-c' line are shown 
on the panel (II). One can see two distinct maxima on 



the profiles. And what is really important is that the 
spatial extent of perturbation is obviously about 80 A. 
The 40 A estimation from Fig. [2] suffers from unsufficient 
contrast of LDOS{x^y)\ev=o image. 

Cross-section of LDOS{x^ eV) along eV coordinate (d- 
d', e-e', f-f and g-g' lines Fig.|4jl)) corresponds to point 
spectroscopy LDOS{eV) dependencies (Fig. QlII)) at 
points of LDOSix) profile, marked by vertical arrows 
in Fig. [4]^II). These are points shown in Fig. [l] by dots 
and arrows. 

Curves d-d' and f-f are taken in between dimers in tt- 
bonded row, while curves e-e' and g-g' are taken on top 
of dimers (Fig. |4]^II) , Fig.jl]). For the whole range of bias 
voltage the values of LDOS collected in between dimers 
are higher than that on top of dimers, except for narrow 
interval in vicinity of Fermi energy, where resides the 
split state. This split state contributes to the increase of 
LDOS on top of dimers in 7r-bonded row. Thus, the pro- 
trusion, consisting of few dimers appear on LDOSix^ y) 
(as well as on STM) image. It follows from Fig. ptill) 
that the contrast of protrusion is higher on (a-a') plane 
than on (b-b') plane, and this indeed can be observed in 
Fig-H 
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FIG. 5. LDOS{x,eV) map along b-b' dimer row in vicinity 
of P atom, located at positions 1 on Ge(lll)-(2xl) surface and 
its cross- sections along denoted lines. Note that impurity can 
appear on LDOS image as protrusion, depression, protrusion 
superimposed on depression etc. 

To illustrate the usefulness of LDOS{x^ eV) map the 
set of cross-sections along spatial coordinate is depicted 
in Fig [5] for P donor atom located at position 1. Pro- 
files are slightly low pass filtered to stress the long 
range features, so they look a bit different comparing to 
Fig. [4]^n). When tunneling bias changes, the LDOS{x) 
profile also changes revealing depdf and protrusions of 
different shape. The profile, corresponding to split state 
energy (and to the presence of protrusion on STM im- 
age) is marked by ellipse. One can easily see that the 
amplitude of protrusion at Fermi level is much less than 
the amplitude of features at other bias voltage (see also 
Fig.glll)). 

Also note, that impurity's LDOS image might have 
elongated hillock like shape at positive bias (empty 
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states). It means the protrusion on the LDOS image is 
not caused by charge density effects (hke charge screen- 
ing). 

To the best of our knowledge this fact was never clearly 
stated. In other words, STM image of Ge(lll)-(2xl) sur- 
face^ (as well as Si(lll)-(2xl) surface^ ^) around surface 
defect is dominated by the split state in vicinity of Fermi 
level, although the amplitude of the effect is relatively 
small. 

To give even more insight into the power of 
LDOS{x^ eV) data representation it is drawn in Fig. [g] as 
quasi-3D surface. Height is given on logarithmic scale to 
increase the image height contrast. The value of LDOS 
is coded both by height and by color with lightning. The 
spatial and energetic positions of specific features of tun- 
neling spectrum can be easily deduced from the figure. 
The split state (zoomed in on the inset) is located at 
Fermi level. It has cigar like spatial shape which directly 
reflects in the shape of protrusion on LDOS{x^ y) image. 

It is also obvious from Fig. |6) that split state really 
fills in the whole width of LDOS{x^ eV) spectra. Thus 
we can not completely exclude the possibility of impurity 
induced electronic features overlap between neighboring 
super-cells of calculation. This overlap can introduce 
some difficult to estimate errors in quantum mechanical 
forces evaluation. This is the main reason, why we have 
increased the size of geometry relaxation surface cell to 
the upper available to us limit. 




Quasi-3D view of LDOS(x, eV) 

FIG. 6. Quasi-3D representation of LDOS(x, eV) along b-b' 
TT-bonded row of Fig.[l] Area, containing split state is zoomed 
in to give clear impression about its spatial structure. Fermi 
level is shown as semitransparent plane. LDOS values are 
given on a logarithmic scale to increase height contrast. The 
values of LDOS are coded both by height and by color with 
lightning. 



Let us now finish the overview of data representation. 
Taking stated above into account one can conclude that, 
given LDOS{x^ eV) it is readily possible to estimate the 
outlook of point spectroscopy curves as well as the shape 
of spatial profiles. That is why the results of electronic 
properties calculations for all 8 possible positions of P 
impurity atom in two surface bi-layers (Fig. [2| are pre- 
sented in Fig. [t] as LDOS{x,eV) maps. 



C. Surface LDOS around P donor impurity 

Let us point out the most important features of calcu- 
lated images. 

As we have discussed earlier, the empty surface states 
band tt* is governing STM image formation for Ge(lll)- 
(2x1) surfac d^ ' ^^ ' in the band gap region. Now the same 
concerns the split state. The STM image is dominated 
by split state in vicinity of Fermi energy. 

The noticeable influence of surface states inside CB 
and VB bands can be inferred from Fig. [7j There are 
LDOS peculiarities near the top of empty surface states 
band tt* and at the edges of filled surface states band tt. 
They are imaged as horizontal bright stripes. 

For all P doping atom positions except position 3, the 
split state is located at the Fermi level. When P impu- 
rity is placed at position 3, the split state can be ob- 
served below Fermi energy Fig.[7|;3). Position 2 IS some- 
what specific. In this case the impurity atom is directly 
breaking 7r-bonded chain, and this strongly influences 
LDOS{x^eV) (Fig. [7|^2)) - at almost all possible bias 
voltage values the impurity LDOS image has two well 
pronounced peaks. 

As with other semiconductors, individual impurity is 
visible in STM experiment when it is located below the 
Ge(lll)-(2xl) surface (Fig [7];5-8)). To the best of our 
knowledge we report this for the first time. Albeit it 
should be rather obvious from simple speculations. The 
crystal lattice is disturbed noticeable far from atomic size 
defect^, and this disturb the perfectness of collective tt- 
bonding in a few dimer rows. There is also another pos- 
sibility. The P is shallow impurity in Ge, its ionization 
energy is 13 meV. Therefore it localization radii should 
be large, in 50 A range. The LDOS, observed by STM 
for Ge(lll)-(2xl) surface must have, in particular, two 
contributions. One quasi-lD, coming from surface recon- 
struction, and another one, coming from ionized donors 
with large localization radii. These two contributions are 
superimposed on each other. The resulting STM image 
would be linear structure, caused by (2x1) reconstruc- 
tion, with wide spots, originating from impurities. These 
spots are poorly visible because of perfect screening by 
TT-bonded electrons, but some influence should exist. The 
calculations for impurity deep below the surface are to be 
done in the future and there is need to re-analyze exper- 
imental data thinking in this direction. Also we should 
mention, that the above speculations must be taken with 
great care, as the common sense might often be mislead- 
ing in surface physics. 

The overall behavior of impurity LDOS images 
strongly differs, depending on the position of donor atom 
in the crystal cell. This is observed even for subsur- 
face defects. The same P donor impurity may looks as 
protrusion, depression, protrusion superimposed on a de- 
pression etc. It depends both on the spatial location of 
impurity and applied bias voltage. 

The atomic orbitals in vicinity of surface defects are 
strongly hybridized. This results in up/downward band 



6 









































AE 


m 

= 


D.OOO 







-X ([Oil]) [nm 



FIG. 7. (1-8) LDOS{x,eV) maps along b-b' dimer row in 
vicinity of P atom, located at different positions in subsurface 
layers of Ge(lll)-(2xl) surface. Numbers denote atom posi- 
tion (see Fig.[2|. Energy difference in eV relative to position 
7 as well as split state position are indicated in the figure. 
Area of split state is zoomed in on every pane. 
(1', 6', 8') LDOS{x, eV) maps along a-a' dimer row for the 
case when two vr-bonded rows are disturbed by impurity atom. 
(4") LDOS{x, eV) map for the case of two donor atoms lo- 
cated at positions 4 and 6. Note split state located below 
Fermi level. Axis directions and images size are indicated 
near the bottom of the figure. 



edges "bending". The insets of Fig. [7| with split state 
areas zoomed in with high contrast, illustrate this. Basi- 
cally, the orbital's hybridization leads to specific spatial 
shape of tunnehng spectra LDOS{x, eV) and, in other 
words, to the appearance of local electronic density spa- 
tial oscillations^^ Let us note, these are not charge 
density oscillations, because they are observed in empty 
states energy range (above Fermi level). Spatial LDOS 
oscillation on Ge(lll)-(2xl) surface were the subject of^^ 
work. 

The energy differences, measured with respect to the 
total energy of a system of 2646 atoms with P donor 
atom at position 7, are shown on every pane of Fig. [7[ 
The difference is not very large. At least, we suppose, it 
does not allow to make any conclusions about the most 
favorable position of donor atom. The difference is large 
for donor positions 8. We do not have any explanation for 
huge energy gain for impurity position 8. At the same 



time this energy difference applies to the huge surface 
slab. Due to slightly different atoms relaxation a few 
electronvolts can easily be acquired by the whole super- 
cell. Also it could be, that the thickness of model slab is 
not sufficient. 

The last row of Fig. [7| will be described later on. 




FIG. 8. (I) LDOS{x, y) maps in vicinity of P atom, located 
at different positions in subsurface layers of Ge(lll)-(2xl) 
surface. Numbers denote atoms position. Maps are given for 
zero bias voltage. Profiled along b-b' line are sketched on 
maps. White lines and arrows are marking the positions of 
maxima in dimer row, nearest to the impurity atom. Crys- 
tallographic directions and images size are indicated on the 
figure. Lines a-a' and b-b' are the same as on Fig. [5] (II) The 
profiles of LDOS{x, y) maps along b-b' line on the same scale. 
Note the shift of profile for position 2. (Ill) High contrast 
LDOS{x,y) map for impurity at position 1, illustrating the 
disturbance of surface electronic structure in a few 7r-bonded 
dimer rows (see text). 

The LDOS (and STM) image of individual impurity 
is dominated by the split state at zero (Fig. [T]) and low 
bias (see above) voltage as illustrated by Fig. [sj where 
zero bias maps of LDOS{x,y)\ev=o foi" different donor 
atom positions are presented together with corresponding 
quasi-3D images. The profiles along b-b' direction (the 
same for Fig. [2] and Fig. |8| are shown on the maps with 
equal scales. 

Three things can be immediately noticed from Fig. [S] 
Firstly, one or two 7r-bonded rows are affected by donor 
impurity. Secondly, one or two local maxima are present 
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on the profile. Thirdly, the distance between maxima 
can be one or two dimers along 7r-bonded row ([Oil] di- 
rection). This is shown in Fig. [s] by thin lines and arrows 
and is summarized in Table. [H 



TABLE I. P donor impurity LDOS image properties 



Atom position 


1 


2 


3 


4 


5 


6 7 


8 


2 rows 


X 










X X 


X 


2 maxima 


X 


X 


X 


X 


X 


X 


X 


Num. dimers 


2 


2 


1 


1 


1 


2 


1 



ment, we have performed calculations for two impurities 
located at different positions in the atomic lattice. First 
impurity was fixed at position 6, while the second was 
sequentially placed at positions from 1 to 4. The results 
for P4-P6 pair is depicted in Fig. [7|^4"). One can see 
that split state for impurity at position 4 was shifted be- 
low Fermi level by adding second impurity to position 
6. Thus we proved, that split state location below Fermi 
energy can be observed for different conditions. The com- 
prehensive analyses of donor pairs is out of scope of cur- 
rent paper. 



Thus, P donor impurity at position 1 is imaged as two 
row feature with two maxima in a row and double dimer 
distance between maxima. 

To be absolutely accurate, not two, but a few rows are 
affected by impurity. The situation is the same as with 
spatial extent of LDOS image protrusion. One should 
either increase the contrast of images (see Fig. [sj^III)), 
either use profiles in analyses (Fig. [sj^II)). Reduction to 
two disturbed rows allows to classify LDOS images of 
impurity located at different positions. 

Having this classification, we can apply it to the test 
case. Si(lll)-(2xl) and Ge(lll)-(2xl) surfaces are sim- 
ilar in many senses. It is possible to perform a simple 
check of our results by comparison with Si(lll)-(2xl) 
surfac^. In general the situation with STM imaging of 
individual impurities is much more simple on Si(lll)- 
(2x1) surface. The empty surface states band tt* and the 
valence band VB are separated by about 0.4 eV gap*^. 
Near Fermi level there are no states, available for tun- 
neling, but empty surface states. That is why the STM 
impurity images on Si(lll)-(2xl) are much easier to clas- 
sify. 

In accordance with Fig. [S] and Table. |l] the conclusions 
of authors can be immediately confirmed. In our nota- 
tions: Fig.28P corresponds to P in position 2, Fig. 21]^ - 
P in position 4, Fig.2b'^- P in position 5. The remaining 
unclear feature (Fig.2d^ ) most probably is the STM im- 
age of P atom, adsorbed on the surface. This statement 
is out of scope of present investigation and will be proved 
in the future publications^^. 

Now we can come back to the lowest row of Fig. [7[ 
Images 1', 6' and 8' correspond to cross-sections of 
LDOS{x, eV) scalar field along a-a' plane (see Fig.[2|, 
i.e. along second (along with b-b') disturbed 7r-bonded 
row. As it can be seen from Fig. [7] the split state is also 
present on these images (see also Fig. [s]). Images 6' and 
8' are similar with the main difference being the image 
contrast. LDOS image for the case of P donor placed at 
position 7 is not shown as LDOS maps along a-a' and 
b-b' dimer rows are almost identical. At the same time 
for other cases the difference between a-a' and b-b' maps 
is rather big. 

The only exceptional case among LDOS{x^ eV) 1-8 im- 
ages is the case of impurity at position 3, when the split 
state goes below Fermi level. To check if this situation 
can be reproduced with slightly different atomic environ- 



D. Local tunneling spectroscopy 

The last problem we would like to discuss is the lo- 
cal spectroscopy LDOS{eV) curves. This article was 
written with experimenter's needs in mind, so we will 
analyze spectroscopy curves as if they were obtained by 
STM. There are two approaches to measure the tunneling 
spectra. One is simple I(V) curve measurement at cer- 
tain surface point. It heavily relies on very high stability 
of mechanical system. Basically this is a case last years. 
Another approach is based on averaging of I(V) curves 
above some surface area. It is less sensitive to different 
noise. 
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FIG. 9. Theoretical tunneling spectra obtained by averag- 
ing over A and B areas (left pane) and by point measure- 
ment (right pane) at a, b, c and d points. Insets depict the 
zoomed in part of LDOS{eV) curves around Fermi energy. 
The middle part sketches the surface area above which the 
spectroscopy was performed. Areas A and B are located above 
a-a' and b-b' 7r-bonded dimer rows (see Fig. [2]), respectively. 
Points b and d are located on top of dimers in a row, while 
points a and c lie in between dimers. The band structure is 
shown on abscissa axis. 



The problem is that two approaches can give different 
from the first sight results. In Fig. [9] we present model 
tunneling spectra for different measurement conditions 
for the case of P atom located at site 1. On the left pane 
there are spectra averaged over A and B areas above tt- 
bonded rows. On the right pane spectroscopy curves at 
points a, b, c and d are depicted. Points b and d are 
located on top of dimers in a row. Points a and c are in 
between dimers. Curves A and B are almost undistin- 
guishable, though the heights of protrusions along a-a' 
and b-b' lines (see also Fig. 2]) on LDOS{x, y) image 
strongly differ. Split state is located close to the Fermi 
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level and this part of spectrum is zoomed in on the inset. 
The difference in averaged spectra is on a few percent 
scale, which apparently is not enough to draw any reli- 
able conclusions. Note nonzero tunneling conductivity at 
Fermi energy. 

As to point spectroscopy, we can easily discriminate 
atomic size features. Dimers at similar positions in dif- 
ferent dimer rows give different tunneling spectra (right 
pane of Fig. [9|. Even relatively small difference above 
elevated features along dimer rows is obvious. 

Looking at the spectra obtained by different methods 
we can conclude, that point spectroscopy does not give 
immediate impression on the band structure, while spec- 
troscopy with averaging does. Averaging of two point 
spectroscopy curves, one on top of dimer and another 
in between dimers, will give curve, similar to averaged 
spectroscopy curve. Also note the vertical scales on both 
panes. Averaging significantly decreases the maximum 
value. 

There are no obvious specific points on numerical tun- 
neling conductivity I(V) dependence (nor on its deriva- 
tive) allowing simple determination of band gap edges 
(see Fig. [TJlII)). I.e. having perfectly defined I(V) and 
knowing LDOS we can not determine the band gap, al- 
though this can be caused by very narrow DFT band 
gap. 

Another conclusion that can drawn from local spec- 
troscopy analyses is that it is almost impossible to iden- 
tify individual impurity on Ge(lll)-(2xl) surface relying 
only on the results of local spectroscopy. As we discussed 
earlier, the maps of L DOS s hould be used together with 
local spectroscopy datsP^I^. 

E. Model's limitations 

Let us specify the strong assumptions used in present 
calculations. Some of them are imposed by very big sim- 
ulation super-cell. 

In particular, we have performed the simulation in 
LDA approximation. It is known it gives non-physically 
small values of band gaps. This can be slightly improved 
by the GGA approximation, but real improvements can 
be achieved only with computationally expensive GW 
many body correctionJ^. At the same time cheap "scis- 
sors" method works quite well^^. 

There is no correction for closed STM feedback loop. 
The values of LDOS are calculated on the plane above 
the surface. 

There is no STM tip density of states in our results. 
In our model we can not account for the surface band 



bending. We simply do not have sufficiently thick model 
slab. Our slab is about 15 A thick, and the depletion 
layer at Ge(lll)-(2xl) surface with n-type of bulk con- 
ductivity is almost 250 A thick. The depletion layer 
strongly affects the picture of tunneling for n-type doped 
Ge samples^. The same concerns the Si(lll)-(2xl) sur- 
face. 

That is why our model STM images do not coincide 
exactly with experimental observations, but neverthe- 
less the correspondence is reasonable. All LDOS{x^ eV) 
maps (except position 3) predicts the presence of protru- 
sion on the STM images at zero (and small) bias voltage, 
which indeed agree with experiments^. We did not find 
any substantial difference when explicitly adding charge 
to the impurity atom. 



IV. CONCLUSIONS 

In conclusion, we have performed the numerical mod- 
eling of Ge(lll)-(2xl) surface electronic properties in 
vicinity of P donor impurity atom located near the sur- 
face. We have found a notable increase of surface LDOS 
around surface dopant near the bottom of empty surface 
states band tt*, which we called split state due to its 
limited spatial extent and energetic position inside the 
band gap. We show, that despite of well established bulk 
donor impurity energy level position at the very bottom 
of conduction band, surface donor impurity on Ge(lll)- 
(2x1) surface might produce energy level below Fermi 
energy, depending on impurity atom local environment. 
It was demonstrated, that impurity, located in subsurface 
atomic layers, is visible in STM experiment on Ge(lll)- 
(2x1) surface. The quasi-lD character of impurity image, 
observed in STM experiments, is confirmed by our com- 
puter simulations with a note that a few 7r-bonded dimer 
rows may be affected by the presence of impurity atom. 

More work is needed to clarify if deep subsurface im- 
purity atoms will be visible for STM on Ge(lll)-(2xl) 
surface and to investigate how the different surface iso- 
mers will affect STM images of individual impurity. 
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